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Astract 

A Procedure is outlined that may be used as a starting point for a perturbative 
treatment of theories with permanent confinement. By using a counter term 
in the Lagrangian that renormalizes the infrared divergence in the Coulomb 
potential, it is achieved that the perturbation expansion at a finite value of the 
strong coupling constant may yield reasonably accurate properties of hadrons, 
and an expression for the string constant as a function of the QCD A param- 
eter. 
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1. Introduction 



In recent work [J], Greensite and Thorn described a simple procedure to produce wave 
functions for states that may describe stringlike features for hadrons in QCD. As an 
Ansatz, a wave function was taken of the form 



N 



ip(xx, x 2 , •■•,x N -x) = AY[(j)(u i ) , (1.1) 

i=i 

where 

Ui = Xi — , Xo, xn fixed. (1.2) 

They then propose to use a variational principle: is chosen such that the total energy 
£, defined by 

£ = T kin + ^ Coulomb , (1.3) 
is minimized. Here, v Coulomh is taken to be simply the Coulomb potential, 

yCoulomb = J2V C (Ui) , V C (u)=-^r. (1.4) 



Subsequently, it is proposed to try "improved Ansatze" , and it is suspected that those 
wave functions that give a string like appearance will carry the lowest amounts of total 
energy £. 

Some aspects of this proposal, however, are less than satisfactory. Certainly, perma- 
nent confinement is not at all a guaranteed property. If we consider the complete set of 
elements of Fock space, we must include all those states in which gluons do not form a 
chain, and therefore it is to be expected that unitarity will not hold if we restrict ourselves 
to chainlike states. A complete set of stringlike states could be found by replacing the 
Coulomb potential by a confining potential, but can one justify such a prescription in 
a theory where the classical limit would not exhibit any resemblance to a confinement 
mechanism? 
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Figure 1: 
The Gluon Chain Model 



1 



In this paper, it is attempted to answer this question. Confinement should be regarded 
as a natural renormalization phenomenon in the infrared region of a theory. The way to 
deal with renormalization counter terms in the infrared is not so different from the more 
familiar ultraviolet renormalization. In Section |2|, we explain how to look upon bare 
Lagrangians, renormalized Lagrangians, and how counter terms may switch position in 
the perturbation expansion. 

In Section the counter term expected for renormalizing the Coulomb potential is 
discussed. At this stage of the procedure, fixing the gauge condition to be the radiation 
gauge is crucial. There are no ghosts then, so unitarity of the procedure is evident, as 
long as we take for granted that the ultraviolet divergences can be handled separately by 
the usual methods (the ultraviolet nature of the theory is assumed not to be affected by 
what happens in the infrared). We then explain how to work with such counter terms 
(Section |). 

Working with non- gauge-invariant, non-Lorentz-invariant expressions will be cumber- 
some in practice, which is why we seek for a more symmetric approach. In Section |5|, it 
is shown that gauge-invariant renormalization effects in the infrared can also easily lead 
to confinement. Again, the renormalization counter terms are essentially local but soft, 
they cannot exist as fundamental interaction terms in a quantized field theory, but they 
must arise as effective interactions due to higher loop contributions. 

Technicalities are worked out in Section |6|, and finally, it is shown that magnetic 
confinement could be handled in a similar fashion, although, of course, the standard 
Higgs mechanism is a more practical way to deal with such theories. 



2. Renormalization 

To illustrate our new approach, we remind the reader of the standard procedures for 
ultraviolet renormalization. In that case, the "bare Lagrangian" , £ bare ; is rewritten as 

C(g, A,---) = £ barc + AC - AC , (2.1) 

where one subsequently interpretes the combination 

C Imarm (g tm , A mn ) = £ baro + AC , (2.2) 

as the finite, "renormalized" Lagrangian. Using a regulator procedure such as, for in- 
stance, a lattice, one may allow that £ bare ; the bare couplings g and the bare fields A 
become ill-defined, or "infinite", in the limit where the regulator vanishes (a j 0), but the 
limit should be arranged in such a way that the renormalized coupling parameters g rcn 
are finite, and the renormalized fields A ren are finite operators. 

AC does contain infinities, but these are precisely constructed so as to cancel all 
infinities that one encounters in calculating the higher loop corrections, order by order in 
perturbation theory. In calculating these higher order corrections, one again uses £ renorm ) 
not £ bare ; because, indeed, AC is again needed to cancel the sub divergences at higher 
orders, as one can verify explicitly. 
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It is important to note now that indeed we could have chosen anything for AC, as 
long as we insist that whatever we borrow to include in the "lowest order" Lagrangian 
£ renorm , we return when a higher loop correction is computed. The procedure pays off 
in particular if we manage this way to keep the higher order corrections, after having 
returned AC together with them, as small as possible. 

What does "as small as possible" mean? This is something one may dispute. The 
renormalization group is based on the observation that when features are computed for 
external momenta (and possible relevant mass terms) at a typical scale /x, then the opti- 
mal choice for AC depends on this scale \i. For small scale transformations, the transition 
from one choice of AC to another may not seem to be very important, but for transi- 
tions from one scale to a very much larger or smaller one, it is. The result is that the 
renormalized coupling strength a s depends on /j, and it may vanish (logarithmically) as 
/j — > oo. 

There appears to be no objection against a wider use of such a procedure. If our 
higher order amplitudes exhibit infrared divergences, besides the ultra-violet ones, we 
could absorb them in AC as well. This time, AC is not expected to affect the coupling 
strengths and the field operators, but a quantity that is ideally suited to be renormalized 
by such terms is the effective Coulomb potential. Thus, after already having dealt with 
the ultraviolet divergences, we add to the "lowest order Lagrangian" £ renorm ; a further 
term AC that affects the Coulomb potential. We just borrow this term from the higher 
order corrections. As soon as these are calculated, we will be obliged to return the loan. 

3. The Renormalized Coulomb potential 

In the sequel, we refer to the usual QCD Lagrangian as £ bare ; even if the ultraviolet 
divergences have already been taken care of. AC will then only describe the new, infrared 
counter terms. Thus, we write 

£bare = -1 F ^ F ^ _ ^ D + + £gaugc-fix + £ghost + AC - AC . (3.1) 

The radiation gauge, 

diA = , (3.2) 

is useful if questions of unitarity are to be kept under control, such as here. As usual, 
Latin indices i, j, ■ ■ ■ are taken to run from 1 to 3, whereas Greek indices fi, u, ■ ■ ■ must 
run from to 3 (or from 1 to 4). The color indices are of no direct concern at this stage, 
since all our expressions at this order are diagonal in the color indices, which is why we 
suppress them, for the time being. In this radiation gauge, the kinetic part of the gauge 
field Lagrangian will be 

£ba re = _i ^ Au) 2 + ^ dvAv f = -|(c^.) 2 + \^A 3 ? + \{d 3 A f . (3.3) 

Here, A$ behaves as an ordinary dynamical variable, except for the constraint ( |3.2|) , 
whereas A Q only emerges in a spacelike derivative term; its timelike derivative, OqAq is 
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absent (it cancelled out in (|3.3|)). Therefore, A generates the instantaneous Coulomb 
interaction between all color charges. 

The Coulomb interaction resulting from the bare Lagrangian is 



^Coulomb ( f ) = -f± . (3.4) 



r 



Our procedure will be particularly important when infinite infrared corrections show up. 
This is exactly what we expect in a theory with confinement. We do not intend to 
"prove" permanent confinement here. Such a 'proof, or at least some justification for 
its assumption, is assumed to be given elsewhere. We take for granted now that the 
renormalized Coulomb potential will be a confining one. Typically, what is needed is an 
'effective' potential that is diagonal in the color indices, and takes up the new form 

n/ 

n r ,r( f ) = - + £M > (3-5) 

where g is an effective string constant. Notive that the signs here are chosen in such a 
way that the Coulomb force keeps a constant sign. The presence of a zero in Eq. (|3.5|) is 
of no significance. 

The familiar Coulomb term, Eq. ( |3.4|) , obeys the Laplace equation 



dV Coulomb (r) = Ana s 5\f) , (3.6) 
whereas the extra term in Eq. ( |3.5| ) obeys the Laplace equation 



B\ 6 \ r \) = p- ■ d 2 ^)=-8n 5 3 (r) . (3.7) 
| r\ \ | r| J 

Therefore, in 3 dimensional Fourier space, the renormalized Coulomb potential is@ 

V-(k) = -^- 8 -^, (3.8) 
k 2 k 4 

so that the Laplace equation for the renormalized potential will have to be constructed 
from 

k 4 

V ren = -A7ia s S 3 (f) (3.9) 



k 2 + 2g/a 

(Here, an admittedly ugly mixed Fourier-coordinate space notation was admitted). Be- 
cause of the extra term 2g/a s , the renormalized Lagrangian should contain the following, 
modified Coulomb term: 

C h - + AC = ... + ldjA U - 20/a ] ) djA , (3.10) 

V k 2 + 2g/a s J 

so that AC must be chosen to be 

AC = -±&A - 2g/as djA . (3.11) 
2 3 k 2 + 2g/a s 3 




Figure 2: The one-loop corrections to the gluon propagator: the gluon, the quark and the ghost term 
(the minus signs are due to the Fermi statistics of the last two; the last diagram, the ghost, actually does 
not contribute in the radiation gauge). 



In coordinate space, this takes the form 

/ A£(x,t)d 3 x = -\ J d 3 x J d 3 f djA {x)G{x-x' )d j A (x l ) , 

where G(r) = — ^- e VV^M . (3.12) 
a s I r\ 

Note that G(x — x') is fairly 'soft', both in the infrared (where it decays exponentially), 
and in the ultraviolet (where it is integrable and considerably weaker than the Dirac delta 
function of the canonical term). On the other hand, this term cancels the canonical term 
completely at k — > . 



4. Infinite infrared renormalization 



In the standard perturbative approach, the one-loop correction to the gluon self-energy is 
infrared divergent. The self-energy correction from the diagrams of the type of Figure ^] 
is 

Cg 2 {k 2 8^-k^k v )\og^ . (4.1) 

It would produce a correction term to the Coulomb potential of the form log(l/fc 2 ) times 
that term. The divere; I — > 00 is the standard type that is to be taken care of by 

ultraviolet renormalization. It is the divergence near \k\ — > 0, which goes as 2Clog(l/|fc| , 
that we now wish to regulate. The counter term produced by Eqs. fl3.11| ) and ( |3.12| ) has 
the form 

- r / (k\ v - k^ky) ■ ft 2 = 2g/a s . (4.2) 
k 2 + fx 2 

The coefficient /x could now be required to be adjusted in such a way that Eq. (|4.2|) 
cancels Eq. ( |4.1| ) "as well as possible" . 

Actually, the first order correction is not quite Eq. (|4.1|) , because the first order La- 
grangian was modified by the AC term. The diagram that has to be calculated actually 
is the one of Figure |^, where the shaded area denotes the effects from the confining po- 
tential, Eq. (|3.5| ). The effects this is expected to have on calculating the higher order 
corrections to the Coulomb potential are sketched in Figure f|. 
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Figure 3: The one- loop corrections to the gluon propagator in the presence of a confining potential. 



Although the counter term ( |3.il| ) is sufficient to produce the above effects at higher 
order, one could also think of replacing it by some gauge- invariant, or nearly gauge- 
invariant expression, such as 



AC = +\J F^(x, t)G(x - x')F llu (x', t) d 3 f 



G ^ = ^L-L e -V^:\r\ . 

a s r 



(4.3) 
(4.4) 



Notice that this is only gauge-invariant in the Abelian case. In non-Abelian theories, G 
must be replaced by a gauge-invariant path-ordered expression. 





Fi gure 4: The one-loop corrections to the gluon potential, a) The standard pure logarithm (in the case 
of massless flavors only), b) Estimate of the modification due to the confinement of the gluons inside the 
loop; the singularity at k — > is exected to disappear, c) The counter term, Eq. (4.2). d) The expected 
total one-loop correction, after subtraction of the counter term. 



On the other hand, in Eq. ( f4.4| ), one could actually allow for more general expressions 
such that the higher order corrections to the Coulomb potential cancel out completely. 
This would not turn all higher order corrections to zero, since we limit ourselves to counter 
terms that are strictly local in time (see Eq. ( |4.3| ). 

Note that, in all cases, the AC that was inserted at zeroth order, is borrowed from 
the correction at higher order (typically the one- loop order), where it is returned, so as 
to render these correction terms as small as possible. There is considerable freedom in 
choosing AC , allowing it to be non-local in space, but we do insist that it stays local in 
time. This will guarantee that the Lagrangian we work with will always produce unitary 
S'-matrices at all stages. It also guarantees that, if the first loop corrections are small, 
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also the higher loop corrections stay small, just like the situation we have in conventional 
renormalization procedures, where this can be proven using dispersion relations. Q 



5. Gauge- invariant counter terms 



The procedure described in the previous section had the disadvantage that gauge-invariance 
is lost. Locality of the counter term in time, and the fact that it is cancelled at higher or- 
der, should minimize the damage done, but still it might be preferable to keep the infrared 
counter terms gauge-invariant from the start. Eq.( |4.3| ) is a step in the right direction, but 
then, the non-local kernel (|4.4j ) should also be made gauge-invariant, which is somewhat 
complicated in practice. It would be easier if we could fabricate the counter terms in such 
a way that the lowest order Lagrangian obtains the form 

£ renorm. = £ bare + A£ = 1^^) . ( 5 .1) 



We have not yet worked out the formal scheme of correcting for this modification at 
higher orders, but it might be doable. The point that will be made in this section is that, 
with certain minor modifications from the canonical theory f(x) = x, which only need 
to be sizable at small values of x, theories can be obtained that provide for rigorously 
confining potentials. 

To demonstrate this, we first rewrite the Lagrangian ( |5.1| ), including a possible 
external source, as 

C(A, 0) = -\Z(<P)F^ U - V(<j>) + J,A^ . (5.2) 

As the field was given no kinetic term here, the Lagrange equations simply demand 
that we extremize with respect to the field, so that equivalence with Lagrangian (|5.1|) is 
easy to establish. In the following, we consider the static situation, and by concentrating 
on one color excitation only, we ignore the non-Abelian terms that otherwise could arise. 

To solve the static case, with fixed sources Jo, the complete equations are easier to 
handle if we first solve the equations for the fields . One gets 

diDi = q = J ; (5.3) 
D = Z{<j>)E , Et= -diA . (5.4) 



By adding to the Lagrangian (|5.2f ) a term 

-^(ZdiAo + Di) 2 , (5.5) 
where Dj is a free variable, we see that what has to be extremized is the Hamiltonian 



!jD 2 
Z(0) 



H = D t A t -C = \ 7777Y + V(<f>) + J A , (5.6) 



which is just — C in the static case, with Eq. (|5.3|) as a constraint. 
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Figure 5: The quark-antiquark vortes in a classical model. Shown are the D field lines (whose divergence 
is controlled by the quark sources), and a sheet crosssecting these lines. Only inside the surface £ 
(indicated as a shaded ellipse) they differ appreciably from zero. 



Let us assume momentarily that the minimum energy configuration (in the presence 
of quark sources) is that of a vortex, as is illustrated in Figure |5|. The D field lines may 
spread over a surface area S . In that case, the predominant field strength D is 

/-) quark /Q 

-> £= = , (5.7) 

£ D K ' 

where Q stands for the (Abelian) charge of a quark. The energy per unit of length is 
£ ^ . ,„„,^n . ( ^ U{D) \ 



,_q^ = min (£[/(£>)) = mm I Q-^- I . (5.8) 

If, indeed, Eq. (|5.8| ) has a non-vanishing minimum, then Fig. |5] correctly represents the 
configuration with optimal energy. Note that, in the Maxwell case, one would have: 

U(D) = \D 2 , (Maxwell), (5.9) 

so then there is no non- vanishing minimum, so that the field lines spread out as usual. 
The minimum exists if, for low values of D , its energy goes as 

U(D) I g\D\ , (5.10) 

which, for simplicity, we take to be the \D\ — > 0. For large D, we expect 

asymptotic freedom to come into effect, so that the Maxwellian behaviour (|5.9|) is resumed. 
For simplicity, also, we take Q = 1, so that g = £ strmg . 

Thus, we see that Eq. (|5.10|) is sufficient to obtain the effect of linear vortex formation. 
We now investigate what this may imply for the function f(x) in Eq.( |5.1| ). 

6. An exercise in Legendre transformations 

First, one has to establish the required relation between Z((p) and V{<p). extremizing 
with respect to implies 



s 



1 Z 1 z 



Figure 6: The relation between Z((j>) and V{4>) producing confinement, and the resulting relation 
between x and Z . 



Regarding U as a function of \D and 0, we have 

dU = i£ 2 di + ^d(±D 2 ) + d^ = I d(±D 2 ) . (6.2) 
In the domain where we require Eq.( |5.10 ) to hold, we therefore have 



(6.3) 



V = - \D 2 d(-) « = ±£ 2 Z. (6.4) 



1 dC/ _ d(gD) _ g 
Z ~ d{\D 2 ) ~ d±D 2 ~ D 

From Eq. ( |6.1|) , one deduces 

1 

At large values of -D, we expect the Maxwellian behaviour Q5.9Q , so that there 

dV 

Z-+1 : ( 6 - 5 ) 

Whether or not V itself tends to a constant there, or to infinity, is not very important. 
All in all, we expect the behaviour sketched in Fig. |6|. 

With this result we now return to Eq. (|5.2|), which we rewrite as 



£ = Zx-V , x = -\F^F tlu . (6.6) 

Extremizing this with respect to <fi gives 

dC „ dV_ 
~dZ 

This allows us to eliminate <fi by writing 



in = - = * = -\F,uF, v = \E 2 . (6.7) 



dC = xdZ + Zdx - dV = ZdX ; Z = — . (6.8) 

dx 

Combining Eq. ( |6.4| ) with ( |6.7| ), we find the relation between x and Z . Inverting this, 
using the relation Z(x), we find C(x) , and we see that at 



x = l£2 = V, ^ = 0. (6.9) 
dx 
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Figure 7: Different possibilities for the function f(x) yielding confinement. The behaviour to the right 
of the extremum (solid curve) must be as indicated. 



Apparently, the function f(x), which in the Maxwell case is just equal to x, must now 
have an extremum at the value (16.91). 



This establishes the function f(x) for all x > \q 2 . How exactly the function continues 
for x smaller than that value is of no direct consequence for the confinement mechanism; of 
course / must approach the classical value x at large \x\ . Figure shows the possibilities. 

Finally, let us derive a curve f(x) that would produce magnetic confinement, i.e. the 
Higgs mechanism. We perform a dual transformation: 

Gfj, u = Z{4>) \ e iiva.pFaP — % F{iv ■ (6.10) 

In the source-free case, the Maxwell equation d^(Z F^ u ) = implies 

d a Gp~i + dpG^ a + = , (6.11) 

so that one may write 

G liV = d ll B v -d v Bp, (6.12) 
and the equations are associated to the Lagrangian 

£(5,0) = ^%^-y(0). (6.13) 

z{4>) 

Electric confinement for the F and A fields implies magnetic confinement for the G, 
and B fields. We see that Z is replaced by \jZ . Plugging in the same equation ( |6.4j) , 
V = \q 2 Z , as Z — > 0, we see a prescribed behaviour when y = \G jlv Gf lv = — x = \B 2 is 
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Figure 8: The function f(x) in the case of magnetic confinement. 



positive: 




(6.14) 
(6.15) 

C -> -gj2y = -q\B\ . (6.16) 



The resulting function f(x) = f(—y) is sketched in Figure |[ 



7. Discussion 

Permanent confinement of electric sources is not necessarily a purely non-perturbative 
quantum phenomenon. One can choose between classical, reasonably looking Lagrangians 
of various kinds that produce this effect. These Lagrangians are non-renormalizable, and 
this means they cannot be used to describe the theory in the far ultraviolet. We must 
assume that the soft, infrared divergent correction terms are due to quantum corrections, 
but they do not have to be non-perturbative quantum corrections. The counter terms 
used in the first chapters of this contribution, are only slightly non-local in space, while 
all are local in time. Arguments that are very similar to renormalization group arguments 
can be applied to conclude that terms of this kind, which become dominating in the far 
infrared, necessarily arise when we scale to longer distances. 
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